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RECONSTRUCTION OF STRATIFIED STEADY WATER WAVES FROM 
PRESSURE READINGS ON THE OCEAN BED 

ROBIN MING CHEN AND SAMUEL WALSH 


Abstract. Consider a two-dimensional stratified solitary wave propagating through a body of 
water that is bounded below by an impermeable ocean bed. In this work, we study how such 
a wave can be reconstructed from data consisting of the wave speed, upstream and downstream 
density profile, and the trace of the pressure on the bed. First, we prove that this data uniquely 
determines the wave, both in the (real) analytic and Sobolev regimes. Second, for waves that 
consist of multiple layers of constant density immiscible fluids, we provide an exact formula 
describing each of the interfaces in terms of the data. Finally, for continuously stratified fluids, 
we detail a reconstruction scheme based on approximation by layer-wise constant density flows. 


1. Introduction 


One of the primary means of measuring waves at sea are pressure transducers — recording 
devices distributed throughout the ocean that read the subsurface pressure. There are many 
advantages to this approach: pressure transducers consume relatively little energy, are inexpensive 
to build and maintain, and provide a wealth of in situ data. As they are deployed submerged, 
moreover, they have a minimal impact on marine life and ocean traffic. 

The use of pressure transducers raises several intriguing mathematical questions: How much 
information can really be gleaned about the flow from just the trace of the pressure on the ocean 
bed? Is it possible to totally reconstruct a passing traveling wave from this data? These issues are 
the subject of the present paper. Specifically, we are concerned with the role density stratification 
may play and whether pressure transducers can be used to track internal waves. 

Let us first discuss the situation for constant density fluids. If one takes the linearized free 
surface Euler equations for a small-amplitude gravity wave, it is simple to show that the surface 
profile can be expressed in terms of the pressure trace; indeed, this mapping is a Fourier multiplier 
whose symbol is referred to as the (linear) pressure transfer function (see, e.g., [6l 118, .251 f27] ). 
More concretely, for a sinusoidal wave with wavenumber k. the connection between the the dy¬ 
namic pressure p and the surface displacement t] can be written as 


( 1 . 1 ) 


cosh (kh m ) 
P = 99 cosh(fci) 


where h m is the height of the pressure meter above the bottom, d is the mean depth, and the 
dynamic pressure p is defined to be the difference between the pressure and the hydrostatic 
pressure, that is, 

P = P -Patm Pg{d h m ). 
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For pressure data obtained from the sea-bed, h m = 0 and hence (11.11) becomes 

pg 

^ cosh (kd) 1 ^ 

which, in the zero depth limit d —>• 0, recovers the hydrostatic approximation 
(1.2) p = pgr]. 

However, it is well known that nonlinear effects play an important role in the surface recon¬ 
struction problem for shallow water waves or for waves in the surf zone (see [3 El [39], for instance). 
Experiments suggest that the predictions of linear models such as D may diverge sharply from 
observation [7i). For waves of large or even moderate amplitude, the discrepancy can be signif¬ 
icant 13 Eg. Another issue with the linear theory is that, due to its reliance on Fourier series 
techniques, it can only treat periodic traveling waves. 

These limitations have motivated a recent push to obtain reconstruction formulas that account 
for nonlinear effects and which apply to solitary waves. Nonlinear nonlocal equations relating the 
trace of the dynamic pressure on the bed to the surface profile of a solitary wave were established 
by several groups of researchers (see [HI m EH). Notably, these papers work directly from the free 
surface Euler equations without further approximation. In P, [TO], the authors obtained implicit 
but exact and tractable relations for periodic traveling waves and performed straightforward nu¬ 
merical procedures to derive the free surface from the pressure at the bed. In PI, Constantin used 
a strong maximum principle argument completely different from the aforementioned approach to 
derive some improved estimates for the wave height of periodic traveling waves. 

The above body of work focuses on irrotational steady waves propagating through constant 
density water. Actual waves in the ocean, however, frequently exhibit a heterogeneous density due 
to salinity or temperature gradients. This density stratification generates vorticity within the flow 
as the effects of gravity and inertia are experienced differently by heavier and lighter fluid particles. 
Moreover, a heterogeneous density distribution can permit large-amplitude waves to propagate 
in the bulk even while the surface remains relatively undisturbed. These so-called internal waves 
have been observed in oceanic areas for many years (cf., e.g. [33] . or the compendium in [23]). 
They are known to have important implications for particle transport pHl , 36;, [37], mixing and 
energy dissipation [211123 EH E3, as well as affecting acoustic propagation E3E3E3. 

The vorticity that naturally accompanies stratification poses a serious mathematical challenge: 
nearly all of the irrotational theory is built on tools (e.g., conformal mappings, nonlocal reformu¬ 
lations via Dirichlet-Neumann operators, variational principles, etc.) that do not have obvious 
analogues in the rotational regime. As a consequence, rigorous results on rotational waves have 
begun to appear in the literature only quite recently. Assuming the vorticity is constant through¬ 
out the fluid domain, which is the simplest model of a wave-current interaction, Ali and Kalisch [2] 
derived a relationship in the long-wave limit that gives a direct map between the surface elevation 
and the pressure at the bottom of the fluid and within the bulk of the fluid domain. Later Visan 
and Oliveras [42] obtained an exact reconstruction formula in this regime without approximation. 

At present, no one to our knowledge has discovered how to adapt these ideas to the case of 
an arbitrary vorticity distribution. Indeed, constant vorticity is a very special (though physically 
significant) case that rarely provides much insight into how to attack the general problem. How¬ 
ever, in [19], D. Henry proved that, if the vorticity and pressure are real analytic, then the trace 
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of the pressure on the bed uniquely determines a solitary rotational wave. The analyticity con¬ 
dition is needed because he employs a Cauchy-Kowalevski argument. In that sense, the surface 
reconstruction problem is well-defined even for rotational waves so long as they are extremely 
regular. 

Unique determinability with stratification. Our first results pertain to the question of 
unique determinability of a stratified steady wave from the trace of the pressure on the bed. 
Mathematically, stratification plays a somewhat similar role in the governing equation for steady 
waves as vorticity. This suggests that Henry’s method can be adapted to analytically stratified 
flows. Following this strategy, we prove that, if two stratified traveling waves share the same wave 
speed, analytic streamline density function, and analytic pressure trace on the bed, then they 
must coincide exactly; see Theorem 13.11 and Corollary 13.II 

However, analyticity is a very strong assumption that is particularly ill-suited to stratified 
waves. Water columns in the ocean are frequently observed to have nearly constant density 
outside of thin transition layers called pycnoclines. Thus, as it passes through a pycnocline, the 
density experiences something close to a jump discontinuity. Indeed, a common practice is to 
view these waves as being a layering of multiple immiscible fluids each with its own constant 
density. In this model, the density distribution is piece-wise analytic, allowing us to iterate the 
Cauchy-Kowalevski approach of Henry and determine all of the free surfaces (cf. Corollary 13.2p . 

The more interesting question is whether one has unique determinability for a continuous but 
not necessarily analytic density (e.g., the existence theory in [33] is formulated for density of class 
C l,a ). For this case, Henry’s argument does not seem to apply and so an entirely new idea is 
needed. In Theorem 13.41 we confirm that unique determinability still holds in the Sobolev space 
setting. Our method is based on reformulating the problem as an elliptic equation on a strip which 
can be treated using the theory of strong unique continuation (cf. Theorem 13.21) . In fact, this 
approach can be used for constant density rotational waves as well, hence Theorem 13.41 generalizes 
Henry’s work to the Sobolev regularity regime (see Remark 13.211 . 

Surface reconstruction. The results above show that surface reconstruction from pressure 
problem is mathematically well-defined for stratified waves, even those with a more realistic 
degree of smoothness. As the second contribution of this paper, we present an explicit scheme for 
performing this reconstruction. In order to explain our ideas, first we must discuss the physical 
setup in more detail. 

We are interested in two-dimensional solitary waves with heterogeneous density g, limiting to 
uniform flows upstream and downstream, and traveling with a constant speed c. Shifting to a 
moving reference frame eliminates time dependence from the system and thus we may assume the 
wave occupies a steady fluid domain 


(1.3) 


= {(x,y) € M 2 : — d < y < r](x)}, 


where the unknown function rj is the free surface profile. We call the fluid continuously stratified 
provided that g € (7(H). We also assume that g > 0 in H and that the fluid is stably stratified, 
i.e., y eA g{-,y) is non-increasing. 
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Figure 1. Example with two fluid layers, Q j and ^ 2 - 


Continuously stratified fluids are difficult to analyze directly since they are rotational. We 
therefore begin by studying the layered model wherein H is partitioned into finitely many immis¬ 
cible fluid regions 

N 

(1.4) Q = (J Hi, f U := {(x,y) € Q : < y < rji(x)} 

i —1 


with the density Qi := constant in each layer. Here we index so that rj o := —d, i)jy := r] and 
Qi lies beneath flj+i, for i = 1,... N — 1; see Figure [1] Note that this is implicitly assuming that 
the free surfaces all have graph geometry. 


The first task is to develop a procedure that allows us to recover a layer-wise constant density 
irrotational solitary waves from its pressure trace on the bed. For this, we can exploit the fact 
that the velocity potential p and standard stream function ip s are harmonic conjugates within 
each layer Hj. Since the density is constant, <p and the pseudo-stream function if := yfolifs 
are harmonic conjugates as well. The pseudo-stream function is better adapted to the study of 
stratified flows than the standard stream function. Then, the conformal change of variables 


(.t, y) H- 


1 1 

— <p(x,y), — ip(x,y) 

c c 


transforms each layer Q t to a fixed horizontal strip (cf. Section \2. 41) . At this point, an adaptation 
of Constantin’s method mi allows us to reconstruct the flow in the first layer up to the first 
interface. However, when we attempt to continue the process and compute the flow in the second 
layer, we are confronted with an unfortunate feature of the conformal coordinates: because the 
velocity potential is in general discontinuous over the interface, the conformal variables behave 
badly as one crosses from one fluid region to the next. This necessitates some additional work, 
but can be overcome with a further reparmeterization of the data on the interface, and by paying 
careful attention to the jump conditions there (cf. Section T4.51) . Ultimately, we are able to obtain 
an explicit reconstruction that determines the flow in every layer, as well as all of the free surfaces, 
directly from the trace of the pressure on the bed. 


Finally, we return to the case of continuously stratified waves using an approximation argument. 
As discussed above, the layered model is a convenient idealization of a continuously stratified fluid. 
In a very recent work, the authors showed that it can in fact be made rigorous for periodic traveling 
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waves in a certain small-amplitude regime (see 0). That is, one can show that as ftlfh 

converges to g in L°°, the corresponding layer-wise smooth waves converge to the continuously 
stratified one. On the other hand, solitary waves in this same regime can be realized as the limit 
of periodic traveling waves taking the period to infinity [Uj. Using a diagonalization argument, 
we are therefore able to reconstruct a continuously stratified solitary wave as a limit of layer-wise 
constant waves considered in the previous paragraph; see the discussion in Section 15.21 

2. Formulation 

In this section, we present several formulations of the steady stratified water wave system, each 
one tailored to one of the problems we later consider. Their equivalence is fairly straightforward 
to show (cf., e.g., [SI Lemma A.2] where it is proved for much weaker regularity solutions). 

2.1. Eulerian formulation. Consider a two-dimensional wave in a stratified water propagating 
to the right with speed c > 0. We assume that the only restorative force is gravity. Adopting a 
moving reference frame, we can view the system as time-independent. The fluid domain U is a 
priori unknown, but, as in (11.31) . we assume that it lies below the graph of some smooth function 
r] and above an impermeable ocean bed at { y = —d}. 

Let us now formulate the governing equations for a multi-layered wave, which includes the 
family of continuously stratified fluid as a special case. With that in mind, suppose that U is 
organized into immiscible strata as in ([1.41) . Then the wave is described mathematically by a 
velocity field (u,v) : U —>• M 2 , a density function g : U —>• M + , and a pressure P : U —>• R. We 
follow the convention that, for a quantity / with domain U, fj : = f |n ; . 

We are interested in (u,v, g) for which each restriction ( Ui,Vi , gi) is reasonably smooth in each 
Uj, but we do not impose any continuity assumptions on U. By contrast, P will be continuous 
throughout the fluid domain. Note that at this stage we do not require that Q{ is a constant. 

To represent a traveling water wave with wave speed c, ( u , v, g, P) must solve the incompressible 
Euler system: conservation of mass 

(2.1a) (u — c)g x + vg y = 0, in Uj, 

conservation of momentum 

(2.1b) g(u-c)u x +gvuy =-P x , in Uj, 

(2.1c) g(u — c)v x + gvvy = —P y — gg, in Uj, 

and incompressibility 

(2.Id) u x + v y = 0, in Uj. 

Here g the gravitational constant of acceleration. We say a fluid is stably stratified if 
(2.2) y i-a g(-,y) is non-increasing. 

The motion of the interfaces is driven by the kinematic boundary condition: 

(2.3a) v = (u- c)d x r)i, on {y = rn(x)}, 
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Note that this includes the air-sea interface. Also implicit above is the fact that 


Vi _ Vj+l 

Hi C c 


on {y = rn(x )}, 


which follows from the continuity of the normal component of the velocity over the interface. 
Since the ocean bed is impermeable, we must have that 

(2.3b) v = 0, on {y = —d}. 


Finally, for the pressure we require 
(2.3c) P = Patm, on {y = rj(x)}, 

where F atm is the (constant) atmospheric pressure. This choice effectively neglects surface tension. 

It will be important for our later reformulations to assume that there is no horizontal stagnation 
in the flow: 

(2.3d) u — c < 0 in fh 

Looking at ()2.lbl) (12. lcll . it is clear that (|2.3dl) prohibits a certain degeneracy in the system. 

In this work, we mostly consider solitary waves. This means that we impose the asymptotic 
conditions 

(2.4) (u,v) —)• (u,0), rj(x) —>• 0 as \x\ —> oo. 

Here u = u{y) is a given function. In most cases, we will assume that u = 0, i.e., there is uniform 
velocity at infinity. For the time being, though, we state things in their full generality. 


2.2. Stream function formulation. The Euler equations are a complicated nonlinear system. 
It is often far more convenient to reformulate them in terms of a scalar quantity. Observe that 
the conservation of mass ()2.1al) and incompressibility (12.ldh ensure that we may define a function 
= ip(x,y) by 

(2.5) ip x = -y/gv, ip y = y/ g(u - c) in Vt. 

This is the pseudo (relative) stream function for the flow, but throughout this paper we will simply 
refer to it as the stream function. From (12.51) and (I2.3dl) , we see that the no stagnation condition 
translates to 

(2.6) ip y < (1 in Q. 

The level sets of ip are called the streamlines. In fact (12.3a|) - (l2.3bl) directly imply that the free 
surface, internal interfaces, and ocean bed are all streamlines. Since (|2.5j) only determines ip up 
to a constant in each 0j, we may take ip to be continuous in H, and set ip = 0 on the air-sea 
interface. Then ip = — po on the bed {y = —d}, where po is the (relative) pseudo-volumetric mass 
flux: 

rn(x) _ 

Po ■■= / Vo(xpy) [u(x, y) - c] dy, 

J-d 


(2.7) 
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which is a (strictly negative) constant. Likewise, we define 

rvi(x) 

(2.8) pi := / 


' —d 


Vq( x iv) i u ( x > y ) - c ] dy, 


so that {y = rji(x)} coincides with the level set {ip = —pi}. 

The conservation of mass (I2.1a,|) implies that there exists a function p : [po, 0] —>• M + such that 

(2.9) e(x,y) = p(-ip(x,y)) in SI. 

We refer to p as the streamline density function and treat it as given. For solitary waves, one 
determines p by examining the flow upstream or downstream. 

Conservation of energy can be expressed via Bernoulli’s theorem, which states that the quantity 


( 2 . 10 ) 


E:=P + | ((u - cf + v 2 ) + ggy, 


is constant along streamlines. This allows us to define the so-called Bernoulli function j3 : 
[0, M ^ M by 


( 2 . 11 ) 


dE 

— (x,y) = -f3(ip(x,y)), 


in Q. 


The dynamics in the interior of each layer is then captured by Yih’s equation 
(2.12a) Aip — gyp'(—ip) + fd{pf) = 0 in Q. 

Note that if in some layer Qj the flow is irrotational and pi is a constant, then (12.12a[l reduces to 
Laplace’s equation for ip. 

Lastly, for the stream function formulation to be equivalent to the full Euler system, it is 
necessary that E satisfies certain jump conditions on the free surfaces. On the air-sea interface, 
P = -Patm and hence (12.1011 gives immediately that 

(2.12b) \Vip\ 2 + 2gg N (y + d) = Q, on {y = r/(x)}, 

where Q = 2{E + gd)\ y=ri . Likewise, for the interior interfaces we use the fact that P is continuous 
in 0 when evaluating E from above and below on {y = jy(x)}. This results in 

(2.12c) \\7-ip i+1 \ 2 - \V-ipi\ 2 + 2g(p i+1 - pi)(r]i +d) = Qi on {y = ru(x)}, 

with Qi := 2{E i+ i + gp i+1 - Ei - gpi)\ y=Vi . 

For solitary waves, the form at infinity strongly constrains the global structure. First, observe 
that (12.41) implies that 

(2.12d) 'V ± ip = yfg(u — c, v) —>• y/p{u — c, 0), rj(x) -A 0 


as x 


oo. 


Let us focus on the simplest case where u = 0. Then, letting h(p) denote the limiting height 
above the bed of the streamline {ip = —p}, from (12.511 we have 


(2.12e) 


o fP 

h( P ) = \ 

Jpo 


'\W 


: ds. 
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In fact, this couples the wave speed and density because we can evaluate the above identity on 
the air-sea interface to get 


(2.12f) 


d = 


f 

J Po 


'vpi} 


ds. 


Now, it follows that 


(2-12g) 


Q = QNC 2 + 2gg N d, Qi = (g i+1 - g t ) (c 2 + 2 gh 


Finally, taking the upstream limit of (12.101) . we can actually compute [3 directly in terms of h: 

1 


(2.12h) 


Pip) = p\p) 


-c 2 + g(h — d) 


This is a well known fact about solitary stratified waves, but an elementary derivation is given in 

ra¬ 
in summary, the stream function formulation of the problem is to find if) continuous in O that 
satisfies Yih’s equation (12.12aj) in each flj, along with Bernoulli boundary condition (12.12bl) on 
the air-sea interface, the jump condition (12.12cl) on the internal interfaces, and the no stagnation 
condition (12.3dl) . In addition, for a solitary wave, we must have that if: exhibits the asymptotic 
behavior ()2.12dp . which in turn implies that the constant Q, Qi are given as in (|2.12gp. 


2.3. Height function formulation. The stream function formulation has the advantage of being 
a scalar equation, but it is still posed in an a priori unknown domain. For that reason, it is useful 
to consider changes of variables that map 11 to a fixed domain. Of course, one should expect this 
to come at the cost of additional nonlinearity in the resulting PDE. A particularly convenient 
choice is to adopt semi-Lagrangian variables 

0 ,y) ( q,P ) := (x, -'*/’)■ 


These coordinates have a long history in the literature on existence of stratified water waves; 
their use dates back to Dubreil-Jacotin m • Note that, because the internal interfaces and free 
surface are level sets of ^, each layer f \ maps to a strip IZi := {{q,p) ■ p G (pi-\,pi)}. Similarly, 
1Z := Ui TZ is the image of 0. 


Let 

h = h(q,p ) :=y + d, 

which denotes the height above the bed of the point with x-coordinate q that lies on the streamline 
{if) = — p}. Notice that h is continuous in 1Z , though its derivatives will in general not be. 

One can show that Yih’s equation (I2.12aj) is equivalent to the quasilinear elliptic problem 

(2.13a) (1 + h 2 q )h pp + h qq h 2 p - 2h q h p h pq - g(h - d)p p hl = -hp/3(-p) in K. 

The Bernoulli condition on the top ()2.12bl) likewise becomes 

(2.13b) 1 + h 2 q + hp{2gph — Q) = 0 on {p = 0}, 

while the jump conditions on the internal interfaces translate to 

1 + (dqhi+i) 2 1 + ( d q hi ) 2 


(2.13c) 


(d p h i+ if 


(d p hi) 


T 2g(pi+i Pi )h — Qi 


on {p = pi}. 
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Here we are adapting our previous convention and writing hi := h \pi i . Finally, simply from the 
definition it follows that 


(2.13d) 


h = 0 on {p = po}. 


The derivation of this system can be found in many places (see, e.g., mm)- 

Restricting our attention to solitary waves with uniform velocity at infinity, (12.1311 can be 
simplified by taking /3 of the form (12.12hp . and Q, Qi as in ( |2.12gp . For convenience, we record 
here the resulting system for a continuously stratified fluid (the layer-wise smooth case being a 
straightforward generalization): 


(2.14) 


(1 + hq)hpp + hqqhp 2 hqhphpq — hpPp 


be 2 + g(h — h ) 


1 + h 2 + h 2 (2gph - Q) = 0, 
h = 0, 


Naturally, the dehnition of h (I2.12cll implies that we also have 
(2.15) h -A h as |g| -A oo. 


in 1Z, 


on {p = 0}, 
on {p = po}- 


2.4. Reformulation using conformal coordinates. Recall that ij} is given by (12.51) along with 
the proscription that i/j is continuous on 0. It follows that ifi can be written 

ip(x,y) = ~Po + ( Vp(x, z ) [u(x, z) - c] dz. 

J-d 

From this expression, it is easy to confirm that if) € C 0,1 (f2). Observe also that Yih’s equation 
(12. 12af) implies that, if (3 takes the form (I2.12hl) . then in any layer where Qi is constant, the 
corresponding velocity field is irrotational and is harmonic. 

We may therefore define a harmonic function (p with domain |J - f2j by taking to be a 
harmonic conjugate of V'lrii- Thus \7p = ^yp(u — c,v), i.e., p is a pseudo velocity potential in each 
layer. This only determines p \up to a constant. Anticipating an imminent change of variables, 
we choose to normalize so that 


(2.16) p(0,—d) = 0 and {(x,y) G H : p(x,y) = 0} is a connected set. 

Note that there exists a unique continuous extension of p\^ to fbut in general p will not be 
continuous over the interior interfaces. This can be viewed as a consequence of the fact that the 
tangential velocity may jump between layers. 

As p and ij} are harmonic conjugates in each fl*, one can employ the hodograph transform: 


X : (x,y) ha ^~p(x,y), ~i/)(x,y) \ =: (t,s). 

Then the image of each fluid layer is a strip 

% := {(t,s) G M x (sj_i,Sj)} 

with Si := pi/c. The fluid domain fl is thus mapped to an IV-slitted strip T := |J i %. 


















10 


R„ M. CHEN AND S. WALSH 


To keep the various formulations clear, designate the height above the bed in the (t, s)- 
coordinates by H: 

H(t,s ) = y + d for (t, s ) € T. 

Therefore the z-th interface {y = r]i(x)} can be parametrized by t: 

(2.17) Y i (t) = H i (t,s i )-d. 

In order to return to the physical variables (x,y), it is also necessary to find a parameterization 
Xi = Xi(t , Si) in a way such that 

(2.18) (. x,rn(x )) = (xi(t,Si),Yi(t)). 


To reformulate the problem in the conformal coordinates, we begin collecting some simple 
change of variable calculations: 


(2.19) < 


dr = 


d y = 


\f0i 


c 

y/Oi 


[(c - u)d t + , 

[~vd t + (c - u)d s ], 


H t =yt = ~ 

H s — Vs — 


cv 


y/ei[(c-u) 2 +V 2 ] 
c(c — u) 


= -X, 


in %. 


From this it also follows that 


( 2 . 20 ) 


d t = H s d x + H t d y , 
d s = -H t 8 x + H s d y , 


v = — 


c — u = 


yfei[(c-u) 2 +v 2 ] 


cH t 

y/Oi {Hf + H 2 ) ’ 
cH„ 


= x t , 


in %. 


Therefore the height function satisfies 

( 2 . 21 ) 


A(t,s)H = 0 , 

H = 0, 


y/0i{H? + H- 2)’ 

in T, 
ons = sq, 


(. H | + H 2 ) (c 2 - 2 gH + 2 gd) = — on s = 0. 

Bn 

For localized solitary waves as considered in [3)0, [H8l PS771 PTTj . we have the following asymptotics 

1 


( 2 . 22 ) 


lim V M H = 0, 

\t \—>-oo 


yfOi 


uniformly for s € (sj_i,Sj). 


lim H(t.,Si) = H(si), i = 1,- ■ ■ N, 
„ |£|—»-oo 


where H(s) = h(cs). Moreover using (12.12e[l to evaluate H at s t yields 


(2.23) 


m Si )= = 

Jpo cvw) y 


Pj +1 ~ Pj _ \ "_ 

^ VbJ^ VbJ^ 
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3. Determinability 


The objective of this section is to establish that the pressure reconstruction problem is math¬ 
ematically well-defined. That is, we prove in various regimes that the wave speed, streamline 
density function, Bernoulli function, and the trace of the pressure on the bed together uniquely 
determine a steady stratified water wave. In particular, among all solitary waves with a given 
form at infinity, to each pressure trace there corresponds precisely one surface profile. 

The main insight underlying the argument is that, due to Bernoulli’s law, the trace of the 
pressure on the bed gives Neumann data for the height function. Since we already know that the 
height is 0 on the bed, together this amounts to Cauchy data. 

To see this, first observe that, in the semi-Lagrangian coordinates, Bernoulli’s law (|2.10l) be¬ 
comes 

1 + h'f 

( 3 . 1 ) E = P + —-l + gp ( h - d )' 

Evaluating this on the flat bed and using the definition of the Bernoulli function P (I2.11|) leads to 

fPO 

(3.2) E\ y= _ d = E\^ =po = E 1^=0 - / P(p) dp. 

Jo 

As before, the form of the flow in the far field (12.411 allows us to compute that the Bernoulli 
constant E\^ =0 is 

(3-3) F|v>=o = E\ y=1) = -Patm + 2 / 3 (^)( c ~~ '^(0)) 2 - 


Therefore, from ()2.13dl) and (j3.1j) - (j3.3p . we infer that 

1 1 f Po 

P(q,Po) + TTHTt -7 - 9P{Po)d = Patm + ~p{ 0)(c - u(0)) 2 - / fi(p) dp. 

2h^{q,p 0 ) 2 J o 

Finally, the absence of stagnation (12.3dl) implies h p > 0, so we conclude 


(3.4) h p (q,p 0 ) = 


2 ( Patm + 


fPo 

-p(0)(c - h(0)) 2 + gp{po)d - / Pip) dp - P(g,po) 

Jo 


- 1/2 


Thus the normal derivative of h on {p = po} is known as soon as the streamline density p, 
Bernoulli function P, limiting downstream/upstream horizontal velocity on the free surface, and 
the pressure data on the sea-bed are specified. 


3.1. Analytic setting. We have seen that, treating the pressure trace as given, ()2.13dl) and (13.41) 
amount to Cauchy data. It is therefore useful to transform the elliptic height equation ()2. 13f) into 
a first-order system that can be approached via Cauchy-Kowalevski. Define 

(3.5) P : = 7"v G: =t-> H:=h. 

rip rip 

This is an abuse of notation in the sense that H above is not the same as the one appearing in 
Section El Nonetheless, for the purposes of the discussion in this subsection, there is little risk 
of confusion. 
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From ()2.3dl) and the identity l/h p = yfip(c — u ) we know that these quantities are well-defined. 
We can then reformulate the height equation (12.1311 as 


(3.6) 


, _ F \GF q + FG q + Gfl(-p) - Ggp p (H - d) 
' p G [ F 2 + G 2 

- Y _ GF q + FGq + G/3(—p) — Ggp p {H — d ) 

t p ~ f 2 + G 2 ’ 



1 

G’ 


with the following “initial” conditions posed at p = po 


Gq 

G ’ 


F(q,Po) = H(q,p 0 ) = 0. 

Here G(-,po) is uniquely determined by the pressure trace P(-,po) via (13.41) . 

Theorem 3.1. Let h E C 2,a (JZ) be a solution of the height equation (I2.14|) representing a soli¬ 
tary water wave with a given stable real analytic streamline density stratification, a real analytic 
Bernoulli function, wave speed c, and downstream/upstream horizontal velocity on the surface. 
Assume that the wave speed exceeds the horizontal fluid velocity u < c through out the fluid. Then 
h is real analytic and is uniquely determined by the pressure function P on the flat bed. 


Proof. The Cauchy-Kowalevski theorem applied to system (|3.6D implies the existence and unique¬ 
ness of F, G, H, which are real analytic in some open neighborhood A f of the bed {p = po}- Hence, 
in view of (|3.5I) . h,h p ,h q are also real analytic in J\f. On the other hand, under the regularity 
hypotheses above, all solutions of the height equation (12.1311 are real analytic in 1Z (cf. [48]; Theo¬ 
rem 5.1]). Therefore the unique continuation of real analytic functions ensures the uniqueness of 
the solution h to (12.141) with Neumann data on the bottom given by the pressure from (13.411 . □ 


In contrast to the case of solitary waves, a fairly robust existence theory for large-amplitude 
periodic steady stratified waves has been developed in recent years (cf. [44, TFJ (2H 2Q, 451146] h 
Since the Cauchy-Kowalevski theory is local, it can be applied just as well on a periodic domain. 
Thus, in order to adapt the argument of Theorem 13.II to periodic traveling waves, it remains only 
to show that the pressure trace on the bed furnishes Neumann data for h there. Unfortunately, 
because this setting lacks an upstream state, we will need some additional information. 

More precisely, observe that by combining (13.111 with (13.211 for a general Bernoulli function fl, 
we obtain the identity 


(3.7) 


h p {q,Po) = 


rP o 


2 E |^ =0 - / flip) dp + gp(po)d - P(q,p 0 ) 


- 1/2 


So long as we know E |^,=o, we can proceed as before. This leads to the following result. 


Corollary 3.1. Let h E C 2,a (7Z) be a periodic solution of the height equation (12.131) for a stable 
streamline density function p E C^ffipo, 0]), Bernoulli function (3 € C^QO, |po|])> and with wave 
speed c. Then h E C U (TZ). Moreover, h is uniquely determined by the trace of the pressure P on 
the bed and the value of E on the free surface. 
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Remark 3.1. There are many situations in which it is possible to find E\^ = q from other infor¬ 
mation about the flow. For instance, any one of the following is enough to determine E |^, =0 : 

(i) the value of u at the crest and the height of the crest; or 

(ii) the value of u at the trough and the height of the trough; or 

(iii) the average value of the kinetic energy (u — c,v) | 2 on the surface. 

For (i), we just evaluate Bernoulli’s law at the crest (where v = 0). Likewise, for (ii), we evaluate 
Bernoulli at the trough. For (iii), we take the mean value on the surface and use the fact that p 
has mean 0. 


These determinability results can be easily extended to layer-wise smooth flows. Note that 
the interior interfaces are streamlines, hence they correspond to straight lines {p = pt} in the 
semi-Lagrangian variables. The first-order system (|3.6I) will hold in each strip TZj. Beginning in 
the strip directly above the bed, we can argue as in Theorem 13. II and Corollary 13.11 to show that 
h G C u (lZi) and is determined there uniquely by the prescribed data, h is continuous over the 
interface and indeed q eA h(q, •) is real analytic. On the other hand, the jump condition (|2.13cp 
allows us to determine d p h 2 (-,pi). Thus we have a full set of Cauchy data for the first-order 
problem in 7^2- Applying Cauchy-Kowalevski once more allows us to conclude h G C^ilZ 2 ). 
Iterating this procedure, we get the following result. 

Corollary 3.2. Let h G C°’ a (lZ) be a solitary wave solution of (12.131) with wave speed c such 
that hi G C 2,a (TZi). Suppose that for each i, G and f3i G C w ([|pj|, |pj_i|]). Then h 

is real analytic in each layer TZi, and it is uniquely determined by the trace of the pressure on the 
bed. 


3.2. Sobolev setting. In this section, we considerably strengthen the above results by relaxing 
the analyticity requirement to a more physically reasonable degree of smoothness. Our argument 
is built off of the following result on the strong unique continuation for a general elliptic equation 
set in the plane. 


Theorem 3.2. PQ Theorem 1] Given a bounded connected open set Ll C K 2 . Let u G W l ' 2 (Ll) be 
a weak solution to the elliptic equation 


where 


Lu := ( a lJ u Xi + a J u) Xj + b l u Xi + cu = 0 in LI, 


(i) L is uniformly elliptic in the sense that there exists a K > 0 such that 

0 < A" _ 1 |C | 2 < < A'|C| 2 , for all ( G K 2 , and for a.e. x G fl; and 

(ii) a*- 7 G L°°(fi), a\ V G L 9 (0), c G L^ 2 (Ll) fo r some q > 2. 


Then the operator L has the strong unique continuation property in the sense that if a solution u 
has a zero of infinite order at a point xq G Ll, i.e., there exists 5 > 0 such that for every integer 


N > 0 


I B r (x 0 ) 


M 2 dx < C]\rr N , 


for all R < 6, 


then u = 0 in LI. 
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With this result in place, we can now state and prove our main theorem on determinability. 

Theorem 3.3. Let r > 4 and h € W 3,r (1Z) be a solution of (12.141) representing a solitary water 
wave with a stable streamline density stratification p € W 2,r (\po, 0]) and a Bernoulli function 
/3 € W 1,r {[ 0, |po|]), such that the wave speed exceeds the horizontal fluid velocity. Then h is 

uniquely determined by the normal derivative h p on the fiat bed. 


Proof. Let h, h E W 3,r (1Z) be two solutions to (12.141) with identical Neumann data on the bed: 

h p (-,po) = h p (-,po) on R. 

The difference u := h — h satisfies an elliptic equation (see also [33] ) 

Jzfu := (a lj u Xi ) Xj + b l u Xi + cu = 0 in 1Z, 

where 


b 1 
b 2 


c 


and 


— i + K, 


.12 


= a 21 = — hph q , 


.22 


= h 2 

,i p'> 


— hppijlq H - hq) ‘2‘h'phpq hphpq H - hqhppi 


•— hqq(hp + hp) 2 hqhpq + / 3 ( p)(hp hphp + hp ) 9Pp(h d(h))(hp + hphp + hp) 
hqhpq + hphqq 

■= -gpph^, 


U = u p = 0, on {p = p 0 }. 

The above boundary condition guarantees that u can be extended to a W 1,r function defined on 
the half-plane {p < 0} that vanishes on {p < po}. Also, since h, h E W 3,r (1Z), from the trace 
theorem and fractional Sobolev embedding theorem (cf., e.g., [ T51 Theorem 8.2]) we know that 


V 2 h, V 2 h € W 1 / 2 ’ r ({p = 0} U {p = po}) C C 0, 2 r ({p = 0} U {p = po})- 


Similarly, we can make sense of p p and on {p = 0} and {p = po}. This way, the functions 
h, h, p and f3 can be extended outside 1Z while maintaining their regularity and ensuring h p > 0 
in the strong sense. Therefore, from the divergence structure of the operator «5f, we see that 
the extended u also satisfies J£u = 0. This, combined with Sobolev embedding, implies that the 
coefficients in 2*f satisfy the the conditions of Theorem 13.21 and hence 2z? has the strong unique 
continuation property. Note that Theorem 13.21 is a local result. Thus we may apply it on any 
finite strip of the fluid domain to show that the restriction of u there vanishes identically. It 
follows that u = 0 on the entire strip, completing the proof of theorem. □ 


As discussed in the previous subsection, h p {-,pfi) is determined by P(-,po) through (|3.4I) . Thus 
Theorem o can be reinterpreted as a statement about unique determinability of h from the 
pressure trace. 

Theorem 3.4. Let r > 2 and h € W 3,r (1Z) be a solution of (12.141) representing a solitary 
water wave with a stable streamline density stratification p € W 2,r (\po, 0]), a Bernoulli function 
/3 € IL 1,r ([po, 0]). Fix the wave speed c, and the downstream/upstream horizontal velocity on the 
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surface. Then h is uniquely determined by the trace of the pressure P on the bed. In particular, 
the pressure function on the flat bed uniquely determines the free surface profile. 

Remark 3.2. The same argument applies to the regime considered by Henry [19], namely con¬ 
stant density flow with a general vorticity distribution. In that case, the height equation is 

{ (1 + llq ) lipp + hqqhp 2hqhphpq — ^pT( P) i PO O P <1 0, 

1 + h 2 q + h^(2gh - Q) = 0, p = 0, 

h = 0 , p = po, 

where 7 is the vorticity function. Clearly 7 plays a nearly identical role as /3 does in (12.141) . and 
so the arguments above go through with only minor alteration. 


4. Reconstruction for layer-wise constant density irrotational waves 

Up to this point, we have proved that, in a number of physical regimes, the pressure trace on 
the bed uniquely determines a traveling wave in a stratified water. The goal of the present section 
is to develop an explicit iteration scheme that reconstructs the internal interfaces of a layer-wise 
constant density and irrotational solitary wave given the pressure data. 


4.1. Surface reconstruction from pressure. In Section 12.41 we introduced the hodograph 
transform and the reformulation of the height equation in (f, s)-coordinates. Recall that the 
density is constant within each layer and we denote the height function as H = H(t,s ) when 
working in the conformal variables. Reading off the Bernoulli law in the bottom layer from 
far-held data yields 


(4.1) 


Phot + 2[{d t H 1 y+ {8^} 


ggid 


-Patm + / gp(p) dp + 
Jpo 



geid , 


where the second term on the right-hand side comes from inverting the streamline coordinates 
from the Eulerian one at infinity. The boundary condition on Hi, (12.191) . and the fact that H s > 0 
then give a formula for d s Hi(-, sq) 


(4.2) 

where 


<9,Hi (t, s 0 ) 


c 

Vft v/c 2 - Mp 


( 4 . 3 ) P(t) := Ao.(zo(*» -If" gp{p) dp . 

Qi f?i Jpo 

Recall that xo is the inverse of x eA <pi{x, —d). 

The surface reconstruction problem can then be stated as follows: From the form of the wave 
at infinity as described by the data H , d , p, and c, and the pressure trace on the bed as described 
by p, find a parameterization for each of the internal interfaces as well as the air-sea interface. 


Remark 4.1. p is a more convenient quantity to work with mathematically, but will not be 
directly observable in most physical scenarios. However, it can be determined uniquely from 
Tbot> the trace of the pressure on the bed in the original variables. This follows from several 
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observations. First, note that the vertical velocity vanishes on the bed, and thus from Bernoulli’s 
law we find that: 

p x (x,-d) = - 

Here we are using the fact that the absence of stagnation implies <p x (0, —d) = y / p 1 (u — c) < 0, 
which selects the negative branch of the square root. Thus xq is characterized as the solution to 
the ODE initial value problem 


-Patm - Pbot(x) + / gp{p) dp + 


Q\cr 


'P o 


1/2 


< X °^ <fx{xo (t)) 

„ *o( 0 ) = 0 . 


Palm - Pbot(xo{t)) + / gp(p) dp + 
Jpo 



- 1/2 


The condition xo(0) = 0 follows from our choice to take </?((),— d) = 0. Once xo is known, p is 
determined uniquely by (14.311 . 


It is important to note that one cannot hope to choose p arbitrarily and have it correspond 
to the pressure trace of an actual traveling wave. We will require a posteriori information about 
H; that is, we need to make use of the fact that it is known that we are observing a solution, 
and that solutions of the system have a certain regularity and decay at infinity. Otherwise, the 
ill-posedness of the problem cannot be avoided. 

Fix an a £ (0,1). We will assume that the wave to be reconstructed lies in the set 

JF = JF{p, H, c, d) := {H £ C' 0 ’“(T) : H — H , H t exponentially localized in t, 

H a solution of (I2.2ip |. 

The existence theory of Turner mum shows that JF is non-empty if the wave speed is super¬ 
critical (cf. Section [5]). In [0], moreover, Amick and Turner prove that it includes waves of large 
amplitude in the case of a two-layer system. 

With this degree of regularity, our procedure holds in the strong sense; with less decay at 
infinity, we obtain only a distributional representation for the wave. Of course, since H £ JF 
is harmonic in each layer, it is C 00 ^) for i = 1,..., N. In fact, the structure of the equation 
and regularity theory for transmission boundary conditions (cf., e.g., [28l Chapter 16]) imply that 
d t H € C'°’“(T). 


4.2. Jump condition in conformal coordinate system. One of the main technical difficul¬ 
ties for multi-layered flows is that the conformal coordinates behave badly on the image of the 
interfaces. Suppose that F £ C'°(D) and let / := F o Then, in general, / 0 C°(T) since 

fi+l(t,Si) + 

In other words, while each JF is a conformal mapping, JF itself is not a homeomorphism. Of 
course, this stems from the fact that ip may have a jump discontinuity over the interface. Recall 
that we chose in (12.1611 to normalize tp so that 1 ({0}) is a connected set. In terms of the 
hodograph transform, this becomes 

(=l#i+i) -1 (0, Si) = (J£i) - 1 (0, Si), for all* = 1,... ,N — 1, 
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which ensures that, at the very least, functions like / above are continuous along the s-axis. 

Given this issue, it is necessary to examine the jump conditions on the internal interfaces very 
carefully. Consider the situation on the i-th interior interface {s = s*}, for some 1 < i < N — 1. 
Let 

(4.5) X i+ := x i+ i(-,Si), Xi_ :=Xi(-,Si), 

where Xi(t,s ) is the restriction of the Eulerian x-coordinate on %. Likewise define 
<L i+ (x) := <p i+ i(x,r]i(x)), $i-(x) := <pi(x,rn(x)). 

This simply says that Xi± and 3>y- are inverse to each other, i.e., 

(4.6) X i± = X iT ($ iT (X i± )), t = $i±(X i± (t)). 

The relation (14.61) allows us to translate from the parameterization of the interface by the 
coordinates in %+\, to that given by the coordinates in define 

(4.7) r i± {t) := $ l± (X iT (t)), 
then 

(4.8) (J^ + i) _1 (t,Sj) = p^) -1 (Ti_(f),Si), (J^+i) _1 (rj + (t),Sj) = (J^) _1 (t, Si). 

We also mention that (12.161) implies 

(4.9) X i+ (0) = Xj_(0) and r*_(0) = 0. 


Suppose that / is a function defined on %U% + i and let F = foJF be its expression in physical 
variables. Then fi + i(t,Si) corresponds to the value of F l+ \ at the Eulerian point (x,r/j(x)) where 
4>j_|_(x) = t. Chasing definitions, we see that this is the same point at which <f>j_(x) = Tj_(t). 
Therefore, the jump in F over the 7-th interface corresponds to the following difference in the 
conformal variables: 


(4.10) lfj i := f i+ i(-,Si) - 

Hence the jump condition (12.12cj) becomes 


(4.11) 


h} + m 




+ 


e(H-H) . = lelc 


on {s = Sj}, i = 1,..., N — 1. 


4.3. Solution of the Cauchy problem. Our procedure begins with the pressure data on the 
ocean bed, then reconstructs the flow in each layer moving from the lowest layer upward. At each 
step, this involves solving a non-standard elliptic problem: Ft is harmonic in each strip %, and 
we will have both the Dirichlet and Neumann data on the lower boundary. This is essentially the 
Cauchy problem for Laplace’s equation (cf. Section [3]) , which is generally ill-posed. In the present 
case, however, we are not concerned with existence theory — we know the wave exists because 
we have measured it — and the asymptotic decay properties (I2.22P will allow us to make sense of 
the solutions as functions. The ill-posedness issue reemerges if the pressure data is noisy, but we 
forestall that discussion until Section 15.41 

The process of solving the Cauchy problem in each strip can be abstracted using nonlocal 
operators. In fact, these operators are nothing but Fourier multipliers whose symbols can be 
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found via elementary separation of variables. Explicitly, for each i = 1 we define an 

operator Ci by 


(4.13) 


(4.12) Ci(4>i,(j)2) := (i^| s=si ,i? s | s=si ), 

where H is the unique solution of the Cauchy problem 

AH = 0, in % 

H = cj>i, H s = cj) 2 , on {s = Sj_i}. 

That is, Ci takes the Cauchy data on the lower boundary of the (i — l)-st layer {s = s*-i}, solves 
the PDE, and then evaluates the Cauchy data for that solution on the upper boundary {s = s^}. 
For now, we do not specify precisely the domain and codomain of Ci, but rather work formally. 

Taking the Fourier transform of (14.131) in the t-direction, denoting by £ the Fourier variable, 
we see immediately that H := TH must take the form 


(4.14) 
Hence 

(4.15) 


#(£, s ) = (j) i(£) cosh [(s - s*_i)£] + 


<M£) 

£ 


sinh [(s - Sj_i)£] 


' ^i(o«>sh[( ai - ai _ 1 )e] + 


<M£) 

£ 


sinh [(sj - Sj_i)£] 


^£<M£) sinh [(s; - i)£] + 0 2 (£) cosh [(s* - Sj_i)£], 


4.4. Reconstruction of the first internal interface. As a model computation, we derive a 
reconstruction formula for the first interface {y = r/ i(x)} given -Pbot('> s o)- From (14.21) we recover 
the normal derivative of H\ on the bed 


(4.16) 


d s Hi(t,s 0 ) 


c 


Note that both sides above are bounded functions of t and hence can be regarded as tempered 
distributions. Performing a Fourier transform in t yields 


d s H 1 ^,s 0 ) = H- 


yfoisjc 2 - 2p 

in the sense of distributions. Applying (]4.151) . we find that 


(4.17) 

(4.18) 


- s) = sinh [b-^ 
3 s Hi(£, si) = cosh [(si - s 0 )£] T ■ 


yfoi\/c 2 - 2p 


V^T-v/c 2 - 2p 


(0 

> (£), se[s 0 ,si], 

W 


From the definitions of x % , Hi, and JgT, we know that the inverse of is given by 

x = x\ (t,s), y = Hi(t,s) — d. 

This together with (14.5jl and (12.171) yields 

(4.19) (x,rji(x)) = (xi(t,si),Hi(t,si) - d) = (Xi_(t), Yi(i)). 
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Note that there is no reason to add a ± to the subscript for Y\ because s eA y(-,s), unlike 
s i-A x(-, s), is continuous. 

From (|2.19D we have the identity d s Hi = d t x\ . Therefore, (14.171) and (14.18P give 

»= smhl(a r 5 ° K 1 J- (—) (0 - ™«t), 

? 1 v^ivc - 2 P J 

-F{*i-K?) = cosh [( Sl - a 0 )?] X \ t r — \ (?), 

{^mV c2 - 2 p J 

since -T{1}(?) = 27 t< 5(?). Using some elementary identities, we can simplify this further to obtain 


«■*» - (W-"))®- ““pH* 01 ta 

« 4*- - i) “> - - Ml " - {A Cfc, - ■)} “' 

Both (14.201) and (14.211) hold in the sense of (tempered) distributions. Ideally, we would like them 
to hold in the strong sense, and to show this we use the idea of Constantin HB. From (I2.19P and 
(14.21) we know that 

(4.22) [c — u(x, — d)] 2 = c 2 — 2p(t). 


On the other hand, we are working with solutions in the set which means that u(x,—d), 
Vu(x, — d) decay exponentially fast as |x| —>• oo. So, by combining (14.221) and (12.201) . we may 
conclude that p and its derivative are square-integrable. In particular, p(t) -A 0 as \t\ —>• oo. 
Moreover, writing 

c _ 1 = _2p_ 

v° 2 ~ 2p(0 \j(? — 2p (t) (c + yjc 2 - 2p JFfj 

we see that that the left-hand side is likewise in L 2 . 

Observe from (12.191) and (14.51) that X\_ = d s H i(-, si), and therefore the terms on the left-hand 
side of (14.201) and (14.21j) are Fourier transforms of square integrable functions. This guarantees 
that (14.201) and (14.211) are indeed equalities for functions. Putting together (14.191) (14.2ip . we 
obtain a parameterization of the first internal interface 


(4.23) X x _{t) = J =+ f t X' 1 ! cosh [(ai - a 0 )?] X 

YQi J -oo I 


(4.24) Y l (t)=X 


VFi 

_i J sinh [(si - s 0 )C] 


1 


X 


I Vc 2 - 2p 


yfol y \J(? — 2p 
- 1 


- 1 


(?) (0 dtf , 




4.5. Reparameterization. Now suppose we have obtained the parametrization of the i-th in¬ 
terface (Xi_(t),Yi(t)). In order to iterate the procedure given above, we will need to compute the 
function r*- that allows us to deal with the singularity in the coordinates on the interface. It is 
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simplest to do this working in the semi-Lagrangian variables. Recall that the jump condition on 
{s = Sj} (which is also {p = Pi}) is expressed as follows 

1 + (d q h i+1 ) 2 1 + ( d q hi) 2 Qi 

~ 2(a r h t+1 y + 2(e p h t r ~ s(£ ’ i+1 '* i+1 - + T = 0 on {p = 

where 

Qi := 2 [(E i+ 1 - Ei ) + g{Q i+ i - Qi)d} . 

Using the fact that h and h q are both continuous across the interface we can solve for d p h t +\ 


(4.25) dphi+i = y/l + ( d q hif 


Qi ‘^9(Qi +1 8i)hi 4- 


1 + (dqhiY 
( dphi) 2 


- 1/2 


on {p = p^. 


Recall that we have the change of variable relations 


1 / hq I 

— = -1p y = ip x , — = 1p x = <Py 

lip lip 

Thus, from (14.2511 we can recover the full gradient V</? at each point on the i-tli interface: 

2-1 — 1/2 n _,_ _ m . , 1 + (dqhi) 2 ] 1/2 


d x <Pi+i(x,rn(x)) = 
dy<p i+ i(x,rn(x)) = 


1 

d p h i+1 (x,pi) 

d q hj+i(x,pi) _ 

d p h i+1 {x,pi) + (d q hi) 2 


— [l + (dqhi) 

dqhi 


Qi ‘^d(Qi+i Qi)hi T 


(dphi) 2 


(x,Pi) 


Qi ■ 1 Qi)hi T 


1 + ( dqhi )- 


(dphif 


1/2 


0,p») 


Let us now compute r,;_. By differentiating (14.711 and using the identity 

(Xj_|_) — ^Pi (Xi +, Tji (4Q_|_ ) ), 

we obtain 

fj_(t) = X*+(l,77-) • Vipi(X i+ ,7n) 


(4.26) 


= X 


1 + (dqhi) 2 


i+ 


dphi 


cX i+ (t) 


(Xi+ , Pi ) (d s Hi(X i+ (t), Si)) 2 ' 


Here a dot denotes a derivative with respect to t. Integrating, and rewriting in terms of the 
parameterization from the lower layer, this becomes 


Ti-(t) = Ti-( 0 ) + 


r*<+W 


Xi+(o) ( d 8 Hi)(T,Si ) 2 


dr = 


r-X+d) 


lx l+ ( o) pQ_(f')) 2 


■ df', 


where we have used the fact that Tj_(0) = 0 by (14.911 . 

Likewise, an ODE characterizing X l+ can be derived by differentiating the second identity in 
m, giving 

1 = <f>'i + (Xi + )Xi + = ( 1 , r/i) ■ V(pi + i(Xi +1 rn)X i+ 

1 + (d q hi) 2 1 1/2 


= X i+ Jl + (d q hi) 2 


Qi 2,g(Qi+ 1 Qi)hi T 


(dphi 


(X i+ ,Pi) 
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Again, initial conditions for this ODE are fixed by (14.91) . and thus can be defined as the 
unique solution of 


X i+ (t) 






Qi — 2g(gi + i — Qi){Yi + d) + 


Qi (y 2 + (x t .) 2 ) 


- 1/2 


5 Si) 


A i+ (0)=X i _(0). 


In this way, we can determine r,_ and X l+ from (Aj_, YJ, X t _). 


4.6. Iteration. In the previous subsections, we saw how to recover the first internal interface, 
and how to handle the jump conditions there. The next step is to iterate this process so that we 
can reconstruct the entire wave. This is done in two stages: first, we obtain formulas that hold 
in the distributional sense. Second, we use the fact that we are working with H £ 34? to deduce 
that they hold in the strong sense as well. 


4.6.1. Iteration in the sense of distributions. For each i > 1, we define the trace operator T,; with 
domain 34 ? by 

(Hi(; Si), (, d s Hi){Si), (d t Hi)(; Si)) £ J%, 

where 

T*(jr) =: = 3% (1) x J% (2) x J^ (3) . 

Simply put, T, evaluates the traces of the solution H and its first derivatives in the region % on 
the interface {s = s,}. The data we are given to begin the reconstruction process is precisely 

3 % := [ 0,- . C ,0 I . 

V VQiy/c*-Mt) ) 

The iteration procedure is as follows. In the {i + l)-st fluid domain %+\, we are solving a 
Laplace equation with Dirichlet-Neumann data on the bottom boundary {s = Sj} obtained from 
previous step: Take the data ( Hi(-,Si ), (d s Hi)(-,Si), (dtHi)(-, Si)) from the solution Hi in %. In 
order to recast the hodograph transformation in we need to access to Tj_ computed in the 
previous subsection. From (14.81) . the continuity of H and dtH across the interface {s = s*} can 
be expressed as 

Si) — Hi{ji-(t), Si), dtHi + i(t, Si) — dfHi (Tj_ (t ), Sj). 

This allows us to compute d s Hi + \{-, Si) from (14.111) : 


- d t Hi(-,Si) 2 . 


d s H i+1 (-,SiY = 

d t Hi(Ti-(-) 

, Si) 2 +d 8 Hi(Ti-(-),Si) 

2 

(4.27) 

c 2 -2g 

Hi (t/—('), Si) - H(Si) 

- 


+- 

c 2 

(6i+1 Qi) 
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We have therefore found a complete set of Cauchy data d s Hj + 1 , H , + \) for the Laplace 

equation (|4.13|) in the next layer % +1 . Solving this, gives (d s H i+1 )(-,s i+ i), (d t H i+1 )(-,s i+1 )). 

In total, this procedure can be written abstractly as a mapping 


1-1+\ i — 1) • • • > N — 1. 

The nonlocal and nonlinear expression for is given by: 

li +1 (u,v,w) = J 7-1 | cosh [(s i+ i - Si)i]u5!T-(C) + smh 
(4.28) 2 f^u^v^w) = J 7 ' 1 |^sinh[(s i+ i - Sj)£] uorT-{f) + cosh[(s i+ i - Sj)£] fi+i(u,v,w)(£)} 
!■[+i(u, v , w) =d t (uo r,_), 


where 

(4.29) 


1/2 


f i+1 (u,v,w)(-) := 


1 


c 2 + 2 g 


tt(r,_(-)) - H(si) 


— w 


v(tU) ) 2 + ^(r i -(-)) 2 


'(f?i+l i?i) 


and the Fourier transforms are taken in the sense of distributions. 


4.6.2. Iteration in the sense of functions. The above scheme can be improved so that the equalities 
hold in the sense of functions. To do this, we utilize the exponential asymptotics (12.221) to conclude 
that 


(4.30) 


Si) — Hi + i(t, Si) y iL(sj), (t, Sj-|_i) y IL(sj_|_i), 

1 as |t| —y oo, 


d s H i+ i(t,s) 


V Qi +1 


for s € [sj,Sj+i], 


with the convergence being exponentially fast. Thus, the differences between the functions and 
their corresponding asymptotics are square integrable. On the other hand, plugging (I4.30P into 
the first equation of (14.28[) . we find that 


7 \ H i+ 1(-, s i+ i) - H(si) - 


^z+1 

V Qi +1 


(0 = cosh [(sj_|_i - «i)C] F {- £(*)} (0 
sinh [(s i+ i - s,)£] 


+ 




3 s ff i+ l(■,*)- 


\/ £*i+l 


( 0 , 


where 5 s iLi+i(-, s*) is determined from (I4.27p . Further simplifying the left-hand side above by 
using the definition of H in (12.23P . we see that 


7 


^H i+ i(t,s i+1 ) - -ff(si+i)j (£) = cosh [(sj+i -Sj)£] 7 | Hi ( Ti -(-), Si ) --H'(sj)} (£) 

sinh [(s i+ i - Sj)£] 


+ 


jNd a ff i+ i(-, ai )- 


V^+i 


( 0 . 


(4.31) 
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in the strong sense. 

Now, to transform back to the original variables, consider 

x = x i+l (t,s), y = s) - d. 

From the relation dfX l+ \ = d s Hi + \ and (14.151) . we know that 

F{d t x i+1 {-,s i+1 )} (£) = £sinh[(s i+ i - Sj)£] J r {^(r i _(-), s*)} (£) 

+ cosh [(s i+ i - s*)£] 7{d s H i+ i(-, s;)} (£). 

Similarly, using the asymptotics (I4.30p . 

1 


(tj Sj+i) — 


exponentially as |t| —>• oo, 


V Qi+1 

and the identity £sinh(a£)J 7 {l}(£) = £sinh(a£)27r<5(£) = 0, we ultimately conclude that 


7 


(4.32) 


d t x i+ i(-,s i+1 ) - -^== j (0 = £sinh[(s i+ i - Sj)£] J 7 {^(rj-Q, «i) - #(s*)} (6 


+ cosh [(sj +1 - s*)£] J 7 d s H i+ 1 (-, - 


V ffi +1 


( 0 - 


Note that Xr i+ ±\_(t) = .Tj + i(t, Sj+i) and y i+ i (t) = Hi + i(t, Sj + i). Hence the (i + l)-st interface 
can be parametrized by 

(4.33) 

X(j +1 )-{t) = ^ f 7~ l sinh [(s i+ i - s*)£] 7 Si) - H(s i+ 1 ) (£) 


+ cosh [(s i+ i - Sj)£] 7 


i (t, Sj) 


yfOi+l 


(0 (O ^ 


(4.34) 


y i+ i(t) = J 7 1 |cosh[(si + i - Sj)^] J 7 Hi(Ti-(t),Si) - H(s i+ i) (£) 

1 


sinh [( S<+ i - Sj)£] ^ 

£ 


(t, Sj) 


(OKt) + ff(si+i)-d. 


V^*+i - 

An argument like that given at the end of Section 14.41 shows that these are equalities in the sense 
of functions. 

In summary, to iterate in the regular function setting, we replace the trace operator F, in 
the previous subsection by another one which evaluates the difference between the trace and its 
asymptotics 

fi(Jtf) = ^ (d s Hi)(-,Si) - J=, (d t Hi)(-,Si)^ =: J%, i = 1,... ,1V. 

Likewise, we normalized the data on the bed so that it has a Fourier transform defined in the 
strong sense: 

/ c 1 X 


:= 0 , 


yfei^/c 2 - 2p(t) y/Ql 


,0 
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A corresponding reconstruction map for the height function X l+ 1 : i = 1,..., N — 1 

can thus be constructed as follows. 


(4.35) 


Z^i = J 7 1 {the right-hand side of (I4.3ip } , where d s Hi + i(-, s*) is given in (|4.27j) . 


2 i+i = d t x i+ i(t,s i+1 ) - 


\/ £?i+i 


= F 1 {the right-hand side of (14.320 } , 


'-(3) _ 


■*i+1 


= &x 


( 1 ) 


tH+ D 


which the corresponding iteration scheme for the parametrized interface reconstruction is 

r t 


(4.36) 


x (i+i)-(t) = 


t 


+ I dt', 


V Qi+1 


i = 1,..., JV — 1. 


y i+ i(t)=^ ) 1 (t) + i7( Si+1 )-d, 


Composition then gives the (i + l)-st layer directly from the data on the ocean floor (which can 
be determined from the bottom pressure readings): 


Ji+i F+\ oX l o ■ ■ ■ oX\ 


. ‘I//? v 

■ <^6 ~t ati+x. 


and the reconstruction of the air-sea interface can thus be obtained at i = N — 1. 


5. Reconstruction for continuous stratification 

5.1. Overview. We now consider the task of reconstructing the free surface of a continuously 
stratified traveling wave from pressure data. The method of the previous section was to adopt 
conformal coordinates where the governing equation is simply Laplace in a strip. As continuously 
stratified waves are rotational, the hodograph transform is no longer available, and the dynamics in 
the interior are too complicated to permit a similar approach. With that in mind, we will instead 
proceed by approximating the continuous wave by a layer-wise constant density and irrotational 
one. 

Before giving the details of this procedure, let us recall some terminology and important results 
from the literature. The existence of small-amplitude supercritical solitary waves in both the 
continuous stratification and layer-wise smooth stratification regimes can be found in Turner 
m [41], where they are built via a limit of periodic solutions with periods increasing to oo. 
Supercritical here refers to the fact that the wave speed is sufficiently large so that the (generalized) 
Froude number for the flow lies below 1. Turner (and all of the contemporaneous existence theory) 
dealt exclusively with steady waves for which the Bernoulli function (3 has the form (]2. 12hh . As 
we have discussed previously, this is the only possibility if one considers solitary waves that have 
uniform velocity upstream and downstream. Because he aimed to construct such waves as limits of 
periodic solutions, Turner naturally looked for periodic waves that also have this type of Bernoulli 
function. It will be important for our later analysis to note that both the periodic and solitary 
waves in Turner’s theory are waves of elevation, that is, their streamlines lie above those for the 
corresponding laminar flow. More precisely, if we define h as in (I2.12cli . then a wave of elevation 
is one whose height function h satisfies 

(5.1) h — h> 0 in \ {p = po}- 
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Our reconstruction method is based upon the fact that small-amplitude periodic waves with (3 
as in (12.12hl) depend continuously on the streamline density function. This was proved recently 
by the authors in [8]. Briefly, the result is the following. 

Theorem 5.1 (Chen-Walsh [8j). Fix a Holder exponent a € (0, 1), a period L sufficiently large, 
a depth d > 0, and a pseudo volumetric mass flux po < 0. Let p* € C l ' a (\po, 0]) with p*(0) = 1 be 
a stable streamline density function and suppose that h * is a corresponding solution of (12.141) for 
a wave speed c* that is supercritical. If /i* is a wave of elevation and sufficiently small-amplitude, 
then the following statements hold. 


(Al) There is a neighborhood U of p* in L°°([po,0]) such that for any p & U with p(0) = 1 that 
is non-increasing and piecewise smooth, there is a solution (h, p , c) to (|2.14l) such that h 
is a strict wave of elevation and c is supercritical. 

(A2) We have the estimates 

(5.2) II ^ _ II Wper (fl) —^11^ P^IlL 00 

and, 


(5.3) 


c c* 


f 

1 

1 

'po 

IVp 

yffl* 


ds = 0{\\p - /o*||l°°)- 


for some constant C\ depending on />*, h *, and L. Here r := 2/(1 — a) is chosen so that 
W 1,r embeds continuously in C 0,a . 

(A4) Moreover, the pressure trace on the bed converges in the following sense. Let a connected 
set I CC [po, 0]\{pi,... ,pn-i} be given with po E I, and assume that p E C 1,a (I). Denote 
by Phot the trace of the pressure on the ocean bed for the traveling wave with density p, 
and let Pbot* he the trace of P * on the bed. Then 


(5.4) 


Pbot — Pi 


bot* | 


/-t 0,o; / 
'--'per V 


< C*2 (lip — P*||l°°([ Po ,0 ]) + \\p - P*llc 1 ’“(7)) 


where C 2 > 0 depends on the length of I, p*, h*, and L. 


The basic structure of the approximation scheme is intuitive, but several nontrivial technical 
complications arise. Take a solitary wave with continuous stratification that is small-ampltiude, 
supercritical, and a wave of elevation. Assume, as in Turner, that this wave can be realized as 
the limit of periodic solutions of increasing period. Each of these periodic waves can in turn be 
approximated by a layer-wise constant density irrotational flow using Theorem 15.11 Finally, the 
multi-layered flows can be reconstructed from the pressure data as detailed in Section EO 


5.2. Layer-wise approximation to continuous stratification. In order to rigorously imple¬ 
ment the approach suggested above, however, we need to more closely examine the convergence 
of periodic waves to solitary waves. Unfortunately, this also requires working in a new coordinate 
system that can be seen as a normalized version of the semi-Lagrangian variables. Specifically, 
we wish to non-dimensionalize and rescale the streamline coordinate p by mapping 

_ I ,1 f p 1 , 

d' cd J 0 ^(s) S ' 


(q,p) ^ (£,C), where f: 













26 


R„ M. CHEN AND S. WALSH 


For each L > 0, the rectangle 1Z L := {(q,p) € (— L, L ) x (po, 0)} is then taken to 

S L := {(e,C): 5 e (-§,§), C e (-1,0)}, 

with each layer 1Z f being sent to 

T T 1 1 

Si ’■= {(C,C) : £ e (--7> - 7 ), C e (Ci-i,Ci)}. where (,■=—. / —j==ds. 

ad cd J 0 yj p( s ) 

We will continue to denote the full strip by 71. Likewise, 7 Zi, S, and Si are each unbounded. 

The height equation (12.1411 is then recast in terms of a new unknown 

(5.5) w = w(£,Q:= - - -C, 

which will satisfy another quasilinear elliptic problem with the rescaled streamline density function 
(5-6) /3(C) := p(p{ 0) 

appearing as a coefficient. Here, w measures the deviation of the streamlines from their height 
at infinity. One advantage of using this quantity is that it decays upstream and downstream and 
thus lies in a Sobolev space set on S. Note that re is a wave of elevation (|5.1I) provided that 

(5.7) w> 0 in 5 \ {C = —1}. 

Following Turner, to w as in ((5.511 . we associate an energy over the set /C C S 

(5.8) F(w,K) := / pdidC 

Jk 1 + 

Note that the absence of stagnation (I2.3dl) requirement translates to the statement that 1+u;^ > 0 
in these variables, hence the above integral is well-defined. Physically, F describes the amount of 
kinetic energy imparted by the deviation from the far field state. 

Let us now more precisely formulate the class of solitary waves for which we can successfully 
prosecute the scheme outlined above. 

Assumption 1. Consider a solitary stratified water wave represented by ( w s ,p,c s ). We require 
that (i) the wave has finite energy 

F(w s ,S) = R 2 < oo, 

that is sufficiently small; (ii) ||«j s ||wi,oo(s) < 6 for some 5 > 0 sufficiently small; (iii) w s is a wave 
of elevation (15.71) : and (iv) w s can be realized as a limit of periodic stratified waves with the same 
energy. This last statement means that there exists a sequence {(tc m ,p, c m )} of continuously 
stratified periodic steady waves with Wm.d^Wm uniformly exponentially localized, 

||V'U) rn ||tyi> oo (5 m ) ^ *3) F(w m ,S ) = R , 

such that there is a subsequence converging to w s in C°’ a (S k ) for all k € N and c m —>• c. 

From Turner eh, we know that waves satisfying Assumption |T] exist. Indeed, much of the 
solitary wave existence theory for stratified waves is built on exactly this type of limiting con¬ 
struction. Nonetheless, one cannot say that every solitary stratified wave falls within this class, 
hence we must make it a hypothesis. 
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Now let (■ w s ,p , c s ) be given as above with {(w m , p, c m )} the corresponding sequence of continu¬ 
ously stratified periodic waves. Then we may apply Theorem 15.II to approximate each ( w m , /3, c m ) 
by a layer-wise constant density irrotational wave (w^, p N ,c^). From (I5.2[) (15.31) we have the 
estimates 

ll^m U;,7l llcper (7T"*j — ^ll^m w m || Wp£('R, m ) — ^lll/ 3 

and 

|c^ -c m | < C\\p N — p\\l°°, 

where r = 2/(1 — a) and Cj = C\{w m , p,m). Note that using (15.31) . we can adjust p N to ensure 
that = Cm, and thus the multi-layer solutions are (w^, p N ,c m ). 


The difference in the energies can then be estimated: 


\F(w%,S m )-F(w m ,S m )\ < [ 

Js r 


< C 2 \\p N — Ps\\l°° + CsllPsllu 00 II w: 


°N 


| Vm 


N |2 
m I 


1 + d^w. 


N 


Ps 


\Vw r 


1 + d^W r 


N 


d£d( 

— Wm]\w 1 ’ 2 (S rrl ) 


< (C*2 + C'iC'3||Ps||l° o ) \\p N — Ps\\l°°, 


where C 2 = C 2 (\\Vw%\\ La °( S m ) , \\Vw^\\ L 2 ^ Sm) ) and C 3 = ^(HVw^|| L <x, (5 m), ||Vw m || L oo( 5 m)). 
From [8, Lemma 3.7, Lemma 3.10] we know that, 

11 11 L°° (<S m ) ) ll^ u 'mllL 2 (5 m ) — CdlVu^mH^gm)) — C(R,5) 

uniformly in N. The exponential localization implies moreover that these constants are uniform 
in m as well. Therefore 


(5.9) 


\F(wZ,S m )-R 2 \ <C\\p N 


p\\l°°- 


Now for any fixed N, following the diagonalization technique as in Turner m, modulo a 
subsequence, we have 

- N " wf, as m —>■ oo, in C°’ a {TZ k ) for all k£ N, 




oo, 


where (w^,p, c(/) is a layer-wise smooth solitary wave carrying an energy 


F(w i :,S)= lim F(wZ,S m ). 

m—too 


The convergence of the right-hand side is a consequence of the uniform (in m) exponential decay 
of and From (| 5.9 f) we know that 

\F(w?,S) -F(w„ 5)| <C\\p N - pHloc. 

Passing to a further subsequence ensures the convergence of c 1 ^ = c m to some c^d So, for any 
e > 0, taking m large enough gives 

ll^s w s\\c 0 ^('R. k ) — II W m w s llc , o,a(7jfc) + ll^m w m\\ c o, a ^k^ + \\W m W s || c o,a (7j fe ) 

< 2e + C\\p N — Ps\\l°°i 

for all k £ N. Finally, letting m -A oo we conclude 
(5.10) ||u; s — w s\\co,a(Tz^) — C\\p ~p\\l°°i 


for all k £ N. 
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Therefore, given the pressure data for , we can apply the reconstruction procedure of the 
previous subsection to recover w s up to 0(\\p N — p\\ ^oo) in c°’« 

5.3. Convergence of the ocean bed parameterization. The argument above shows that we 
can approximate the continuously stratified solitary wave via a sequence of layer-wise constant 
density irrotational solitary waves. However, its main weakness is that it requires the pressure 
data not for the continuously stratified solitary wave (which is observed), but for a limiting 
(sub)sequence which are only known to exist mathematically. This is partially ameliorated by (15.411 
which guarantees that the pressure traces of the subsequence converge to that of the continuous 
solitary wave. 

But there are two issues remaining that must be addressed. The first of these is rather benign: 
in Section 01 we begin the reconstruction process with the data p. It is therefore necessary to 
confirm that the convergence of the pressure likewise extends to this quantity. Note that because 
p is defined in conformal variables — which do not have an obvious meaning for continuously 
stratified flows — it is not immediately clear that this should be true. The second more serious 
issue is the sensitivity of the reconstruction scheme to noise in the pressure data. This we discuss 
in the next subsection. 

Let h be the height function in semi-Lagrangian coordinates for a steady wave localized near 
the crest (it can be solitary or periodic). Suppose that, as above, there is a sequence {h N } of 
layer-wise constant density steady wave localized near the crest that converges to h. Then for 
each h N , the conformal coordinates are well-defined and we have 

t = -\p N (x% ( t),-d ), 

and hence 

-c = ip” (i), -d)x$ ( t ) = -ipy{x$ (t),-d)x$ (i). 

But then, using the change of variables formulas, this implies that Xq solves the ODE 

Xq = hp( x o,Po), Xq (0) = 0. 

In view of this, we define xo = xo(t) to be the solution of 
(5.11) x 0 (t) = h p (x 0 (t),po), x(0) = 0, 

where the dot denotes differentiation with respect to t. Formally, the initial value problem (15.111) 
is the limit of the ODE satisfied by as N -A 00 . We hope to show that x^ likewise converges 
to x$, which then implies p^ approaches p. 

Proposition 5.1. Let {p N } C L°°([po>0]) be a sequence of piecewise constant stable streamline 
density functions, with p N -a p in L°°, where p € C 1,a (\po, 0]) is a stable smooth streamline 
density function that is constant on some neighborhood I of the bed. Let h and h N be the height 
functions for a traveling wave localized near the crest with period L € (0, 00 ], the speed c and c N , 
and streamline density function p and p N respectively. Assume that 

<V P(Po) = c N ^/p N (p 0 ), 

and that there exist constants M, 7 > 0 independent of L and N such that 
\h P (q,Po) ~ h P (Po)\, \hp (q,p 0 ) ~ h p (p 0 )\ < Me ~ 7|g| . 
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If L = oo, i.e., the waves are solitary, also suppose that h p (-,po) —> h p (-,po ) in L°°(R). Then it 
follows that 

Xq -A xq in L°°(M). 


Proof. Denote z N := xo — x^, so that z N satisfies 

(5.12) z N (t) = f [hp(x 0 (T),p 0 ) - h^(x$(T),p 0 )\dT. 

■Jo 

The no stagnation condition implies that there exists a > 0 such that 

h p (-,po), hp (-,p 0 ) > a > 0, for all N > 1. 
Recalling (15.111) . this guarantees that 

xo(t), Xq (t) > at, for all t € [0,oo). 


Fix e > 0, and choose T > 0 sufficiently large so that 

pOO 


L 


e~ ait dt < —-j—:. 
j' 4 M 


Starting from (|5.12p , we write 

pT poo 

II^IU- < / \hp(x 0 (T),p 0 ) - (x$ (r), p 0 )\dr + / [/ip(x 0 (r),po) - /i^(x^(r),p 0 )] dr 

io Jt 

= : I + II. 

To estimate II, we exploit the fact that h p (po) = h p (po), and the exponential localization: 


II < 


poo 

\hp(x 0 (T),p 0 ) - h p (po)\dr + J \h p (xq (r),p 0 ) - h p (p 0 )\dr 

s: 


<2 Ml e~ iat dt.<-. 


On the other hand, for I we use the simple estimate 

I < T\\h p (-,po) - h p (■,po)||l° o < T\\h p (-,po) - h p (-,Po)||l°° < 
for N sufficiently large. Observe that in the periodic case L < oo, we have 

IIM-.Po) - K (uPo)|| l°° < C (Hp^ - p\\l°°([po,o]) + II P N ~ Pllcu“(7)) > 

and hence the L°° convergence of h p on the bed is implied by our previous assumptions. 

The argument above gives the convergence of Xq to xq in L°°([0, oo)); the argument for 
L°°((—oo, 0])) is similar. □ 
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5.4. Pressure noise. Let us now discuss the sensitivity of the above considerations to pressure 
data. First, it is important to note that any scheme that is attempting to reconstruct from pressure 
must come to terms with the physical limitations of pressure transducers. It is of course impossible 
to perfectly read the pressure, and in particular, existing instruments fail to capture very high 
frequency components. Likewise, computers are finite precision machines, and so any practical 
implementation of the recovery scheme will necessitate additional approximation. For instance, 
it is common to use a windowed transfer functions that simply ignores the high frequency range 
above a certain cut-off. Mathematically, this procedure corresponds to applying a Littlewood- 
Payley projection onto a Unite frequency ball. However, because we are working with an ill-posed 
PDE system, it is very challenging to prove that convergence occurs as the cut-off radius is 
increased. 

This seems endemic to all reconstruction from pressure procedures, even for constant density 
irrotational flows. Nonetheless, actual numerical experiments suggests that the schemes converge, 
and give reasonable agreement with observation mm- 

Unfortunately, the reconstruction process for continuously stratified fluids given above possesses 
an additional source of error. Observe that in our approximation we take the pressure reading 
of a solitary wave with continuous stratification and treat this as pressure data for a layer-wise 
constant density irrotational wave. But, the unique determinability result indicates that the 
pressure trace for these two waves cannot possibly coincide. We can therefore think of this as 
attempting to run the reconstruction scheme of Section |4] with noisy pressure data. 

To be more concrete, denote 

(5.13) Phot = Phot* + Pern 


where -Pbot* is the bottom pressure reading of the solitary wave with continuous stratification, 
and P^ ot is the trace of the pressure for the IV-layered solitary wave. From (14.31) we then have 

(5.14) p N = p + — + — [ g [p{p) - p N (p)\ dp=:p + p err . 

61 61 Jp 0 

In order to recover the true IV-layered flow, one should use p N (t) in the iteration scheme (14.361) . 
rather than p(t). Hence the error from our iteration scheme consists of two parts: one from the 
limiting process by taking large number of layers; the other from the pressure error. 

Let 

From (14.24j) the error in the reconstruction of the first interface is 

(5.15) Y* - Y, = F- 1 j sm M0h-*o)^ [/(p w } _ /(p) ] (0 j . 


The difference term can be computed as 
f(p N ) - /(P) = 


2c 


\JQi (^ 2 - 2p)(c 2 - 2p JV ) (Vc 2 - 2p + \J r 2 - 2p x ) 


• (P“P 


N\ 
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To justify the inverse Fourier transform, or more precisely, to make sense of and Y\ as 
functions, one needs to ensure that T [/(|T V )] and fF[f(p)] have sufficient decay, e.g. 

e (si-so)|Cljr [/(^)] (£), e (^i—o)kl^- [/(p)] (£) € L 2 (R). 

In fact, this can be proved rigorously if the continuously stratified wave has constant density 
in S\. Following Clamond-Constantin [10], we see that the pressure transfer function p N (t) can 
be holomorphically extended to a function p jV (z) with z = t + is defined in the strip {(f,s) € 
R x (2so — si — e, si + e)} for some e > 0. Likewise, the analyticity of p^ (z) and the no stagnation 
condition imply that the function 


/(P A ») = 


- 1 


is analytic in 1 x [—A,/i], where X = fi = s\ — sq + 
e.g. i2, Theorem 4]) implies that 


Vei y y/c 2 — 2p N (z) 

Thus, the Paley-Weinder Theorem (cf., 


e (si-s 0 )|f|jr (£) € L 2 


However the above argument does not extend beyond the region of constant density, and hence 
the decay of ^[/(p)] (£) becomes very hard to verify for more general stratification. Worse still, 
because one loses analyticity crossing the internal interfaces, it is unclear how this reasoning could 
be used in any layer but the first. For this reason, we cannot currently prove that our method 
represents a convergent numerical scheme. On the other hand, as it does succeed an analytical 
approximation method, we anticipate that it will have applications for studying the qualitative 
behavior features of continuously stratified solitary waves. 
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